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Abstract. - Based on the model of F/S/F trilayer with atomic thickness we discuss 
the relative roles of pair-breaking and proximity effects, as a function of the exchange field, 
of disorder and of a finite thickness in the superconducting layer. The exchange field can 
be small (weak ferromagnets) or large (strong ferromagnets) compared to the superconducting 
gap. With weak ferromagnets we show the existence of a reentrant superconducting gap for the 
F/S/F trilayer with atomic thickness in the parallel alignment (equivalent to the F/S bilayer). 
Qualitatively small disorder is equivalent to reducing the value of the hopping parameters. 
In the presence of a finite thickness in the superconducting layer the superconducting gap in 
the antiparallel alignment is larger than in the parallel alignment, meaning that pair breaking 
dominates over the proximity effect. 



Introduction. - Many interesting phenomena take place in mesoscopic devices when 
superconductors are connected to ferromagnets. For example the superconducting pair am- 
plitude induced in a ferromagnet oscillates in space [21|31^i5j and these oscillations can give 
rise to the 7r-couphng. The 7r-couphng manifests itself in superconductor / ferromagnet (S/F) 
multilayers as oscillations of the superconducting critical temperature with the thickness of 
the ferromagnetic layers S/F/S 7r-junctions have been recently probed with various 

experimental techniques P HTUIITT] . 

The proximity effect in F/S/F trilayers has recently focused a renewed interest. It was 
established long ago that with insulating ferromagnets the superconducting transition tem- 
perature is larger in the antiparallel alignment because of the exchange field induced in the S 
layer ^21 that tends to dissociate Cooper pairs. Following this theoretical prediction two ex- 
periments were performed, one with metallic ferromagnets |13| and the other with insulating 
ferromagnets 14 and it was shown in both cases that the critical temperature was larger in 
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the antiparallel spin orientation. On the other hand with metaUic ferromagnets the theoret- 
ical description should also include the proximity effect, namely the possibility that Cooper 
pairs from the superconductor can delocalize in the ferromagnetic electrodes. The proximity 
effect in F/S/F trilayers is unusual because it may involve spatially separated superconduct- 
ing correlations if the S layer is thin enough |15llT]. The basic physics can be understood 
from considering a model of half-metal ferromagnets: in the parallel spin orientation Cooper 
pairs from the superconductor cannot be transfered in the ferromagnetic electrodes. As a 
consequence the zero temperature superconducting gap is not affected by the coupling to the 
ferromagnetic electrodes T . In the antiparallel spin orientation Cooper pairs from the super- 
conductor can delocalize in the ferromagnetic electrodes: the spin-up electron can tunnel in 
the spin- up electrode and the spin-down electron can tunnel in the spin-down electrode 
The superconducting gap is reduced by this proximity effect associated to spatially separated 
superconducting correlations jl5| and the zero-temperature superconducting gap is larger in 
the parallel alignment (Ap > Aap) |lll5j . On the other hand it was shown that the supercon- 
ducting transition temperature with metallic ferromagnets is, like with insulators, larger in 
the antiparallel alignment (TTj {T^ < T^^) which has been probed in recent experiments |18j . 
In a recent work Buzdin and Daumens pp proposed to reconcile these apparently contradic- 
tory results for the gap at T = and for Tc, by calculating the critical temperature of a 
F/S/F trilayer with atomic thickness within the Stoner model, and the zero temperature gap 
within a model of half-metal ferromagnet. They indeed found that < T^^ but Ap > Aap- 
This suggests that the proximity effect plays a dominant role in the determination of the zero 
temperature superconducting gap and that pair breaking effects play a dominant role in the 
determination of the critical temperature. The goal of this Letter is to determine whether 
this picture is robust, including realistic ingredients such as a strong or weak exchange field 
in the ferromagnets, disorder or a finite thickness in the superconductor. 

The model. - Let us start with the F/S/F trilayer with atomic thickness Q]. The 
superconducting layer is described by the BCS Hamiltonian 

WbCS = XI (4,aC/3,<7 + 4,<TCa,<T) + T^.i + KCa,lCa,]^ , (1) 

(a,/3),(T a 

where a and P correspond to neighboring sites on a square lattice. The ferromagnetic elec- 
trodes are described by the Stoner model 

(q,/3),(T i 

where hex is the exchange field. The ferromagnetic layers a and b are connected to the 
superconducting layer by the Hamiltonian 

= J2 (ca,a,a(f,)Ca,cT,S + 4, a.sCa,^.a(f,)) • (3) 

a,cr 

Strickly speaking the individual layers of the F/S/F trilayer with atomic thickness might be 
unstable against charge or spin density wave ordering. We view the F/S/F trilayer with atomic 
thickness as a toy model for more realistic models involving a finite thickness in the F and S 
layers so that we can safely neglect these instabilities. 

The Green's functions of the superconductor are given by ^U] 
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[uj — e{p) + i5] [uj + e{p) — iS] ' 



(5) 



and similar expressions are obtained for g1;\ and /^'J,. The labels "1" and "2" in and 
((SJ) correspond to the Nambu indexes. e(p) corresponds to the quasiparticle energy e{jp) = 
\/aF+^ , with ^ = h^p^ /{2m) — ep the kinetic energy determined with respect to the Fermi 
level. The coherence factors are given by = (1 + £,/e{p))/2 and Vp — {1 — £,/e{p))/2. The 
"11" Green's function of a spin- up ferromagnetic electrode is given by g^'l ~ l/[tj — ^ — /icx + 
iS sgn(^ -I- /lex)] , and a similar expression is obtained for the "22" component. 

Perturbative expansions. - The fully dressed Green's function corresponding to the 
F/S/F trilayer with atomic thickness is determined through the Dyson equation Ga,a = 

ga,a + ga,aia,aga,aia,aGa^a + ga^aia.bQbdbMG a,a- The self-consistcucy relation takcs the 
form JU| 

To order and for strong ferromagnets (Aq <^ h^^) we obtain: 



AoJ 2 



(7) 



where Aq is the superconducting gap of the isolated superconductor, 
ferromagnets (t <i; /lox < Aq) we obtain 



In the case of weak 
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We deduce that in both cases the superconducting gap at order is reduced by the proximity 
with the ferromagnetic layers, independently of the relative spin polarizations of the F layers. 
In the case of strong ferromagnets a perturbation theory to fourth order in t leads to 
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where we introduced a small cut-off rj to regularize the logarithmic divergence of the inte- 
gral over ^. The (positive) difference between the superconducting gap in the parallel and 
antiparallel alignments is of order (i//icx)* and is thus very small compared to Aq. 
For weak ferromagnets a perturbation theory to order leads to 
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In this case there are two small parameters in perturbation theory: I^/Aq and /i^^/Aq. Per- 
turbation theory is well suited for understanding the small parameters in the problem but is 
not quantitatively reliable because of logarithmic divergences in the limit i] ^ 0. 

Numerical simulations of the F/S/F trilayer with atomic thickness. - Non-perturbative 
solutions valid for arbitrary interface transparencies and for finite temperatures can be imple- 
mented numerically. For a F/S/F trilayer with atomic thickness we can solve explicitly the 
Dyson equation and obtain the expression of the dressed propagator G^;^^: Gj^'^ — fa ^/T), 
with 
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Fig. 1 - Variation of the superconducting gap A(T) as a function of temperature T for the F/S/F 
trilayer with atomic thickness (a) for strong ferromagnets and (b) for weak ferromagnets and with 
the parameters h = 1, m — 2, kp = 1, \ = 0.32. The Debye frequency is set to ujd = h^kp/{2m). For 
both figures we use t — ta ~ tb- For strong ferromagnets we used t — 0.01. For weak ferromagnets we 
used h = 0.003 and t = 0.00235 (1), t = 0.00230 (2), t = 0.00225 (3) and t = 0.00175 (4). The open 
symbols correspond to the parallel alignment and the filled symbols correspond to the antiparallel 
alignment. 

where the self-energy is diagonal in Nambu space: Sj^^'^^ = {ta9l]a^ +tb9b^b^^)^Ti,r2J where 
T = 1 , 2 is a Nambu index. 

The temperature dependence of the superconducting gap in the case of strong ferromagnets 
is shown on Fig. n}(a). We obtain Ap(0) ~ Aap(O) but T^^ > . Because electrons in 
the S layer can make excursions in the F layers an exchange field of order (i^ + i^)//iox is 
induced in the superconducting layer in the parallel alignment |12[I20| . This pair-breaking 
mechanism is absent in the antiparallel case. The order of magnitude of the exchange field is 
obtained by evaluating the density of spin-up and spin-down electrons in the superconductor. 
Superconductivity breaks down at zero temperature if the exchange field in the F layer becomes 
smaller than a critical exchange field hi^ of order hi^ ~ (t^ -I- i^)/Ao(0). For the parameters 
corresponding to Fig. ^we obtain from the numerical simulation — 0.092 ± 0.002 which 
is of order 2t^/Ao. The F/S/F trilayer in the parallel alignment with the parameters on 
Fig. nKa) is close to the breakdown of superconductivity and this is why we obtain large 
values of - T^^ . This regime may be of interest from the point of view of realizing a 
superconducting spin valve. Yet, this would require a fine tuning of some parameters such as 
the interface transparencies. 

With weak ferromagnets (see Fig. n-(b)) the most interesting case corresponds to /lox 
slightly smaller than Ao(0). We obtain Ap(0) > Aap(O) for small values oi t = ta ^ h and a 
crossing between Ap(T) and Aap(T). For larger values of t we obtain Ap(T) < Aap(T) for all 
values of the temperature. If the interface transparencies are sufficiently large the gap Ap(T) 
in the parallel alignment is minimal at T = T* and reenters at a lower temperature. This 
behavior is apparently due to a cross-over between two regimes: Ao(T) > /lex for x < T* 
(weak ferromagnets) where the small parameter in perturbation theory is hc-x/ ^o{T) and 
Ao(T) < hex for T > T* (strong ferromagnets), where the small parameter in perturbation 
theory is Ao(r)//iox- Given the perturbative expansion like © we see that the pair-breaking 
effect of magnetism onto superconductivity is maximal for hex — Ao(T). We also verified 
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Fig. 2 - (a) Variation of the critical temperature as a function of disorder in the S layer for strong 
ferromagnets. (b) Variation of the superconducting gap as a function of disorder in the S layer for 
weak ferromagnets. Both cases correspond to the F/S/F trilayer with atomic thickness. We used the 
parameters h = 1, m = 2, kp ~ 1, \ = 0.32. The Debye frequency is set to ujd ~ Ti'^k%/{2m). For 
strong ferromagnets we used t — ta = tb = 0.01. For weak ferromagnets we used h — 0.003 and the 
temperature is set to T = 0.00025. In (a) the open symbols correspond to the parallel alignment and 
the filled symbols correspond to the antiparallel alignment. In (b) the open symbols correspond to 
the parallel alignment. The solid line corresponds to the antiparallel alignment with t — 0.0220 and 
the dashed line corresponds to the antiparallel alignment with t = 0.0235. Disorder is measured by 
for oj > A. 



that for smaller values of hex we obtain larger values of T* which is in agreement with this 
scenario. 

Effect of disorder. - To include disorder in the superconducting and ferromagnetic layers 
we cannot use Usadel equations since we want to describe spatially separated superconducting 
correlations. Instead we use a perturbative treatment of disorder based on Ref. (19||. It is well 
known that the superconducting gap of an isolated superconductor is not affected by non 
magnetic impurities J19I and we could verify this in our simulations. We show here that this 
is not the case in the F/S/F trilayer and that introducing non magnetic impurities has the same 
effect as reducing the value of the hopping parameters. The self-energy associated to disorder 
in the superconductor takes the form E^-^ = nau'^gl;^^{0) , where ria denotes the concentration 
of impurities, Ua is the scattering potential, supposed to be isotropic, and (7^'a(0) is equal to 
the propagator of the superconductor evaluated at i? = 0. The '11' component of the total 

self-energy is given by Ei^i-'"* = + tlgl^^l + tlg^ -f tl (g^)' + tl (gli)\\\ and 

a similar expressions are obtained for 5]^'^''°*. The '12' component of the total self-energy is 
given by I]i'2:tot _ ^^1,2^ ^fiih E^'^ = «ctU^g^'_^(0). E^ and E;, correspond to the self-energies 
due to disorder in the ferromagnetic electrodes. The expression of the '12' component of the 
averaged dressed propagator is given by 

Si'^}, (12) 

where 2? is given by (|ll|l with Et being replaced by Eji^"' . We implement the simplest approxi- 
mation: we neglect the corrections to the tunnel vertex in which a quasiparticle from S scatters 
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Fig. 3 - (a) Variation of the superconducting gap as a function of temperature for strong ferromagnets 
for two values of the thickness of the S layer, (b) Variation of Ap (T')/Aap (T) as a function of the 
hopping amplitude t — ta ~ ti, for weak ferromagnets. We used the parameters h = 1, m = 2, kp = 1, 
X = 0.32, t' = 0.1. The Debye frequency is set to ujd = h^kp/(2m). For strong ferromagnets we used 
t = ta = tb ~ 0.02 and /icx ~ 0.1. For weak ferromagnets we used /lex ~ 0.003 and the temperature 
is set to r = 0.0004. In (a) the open symbols correspond to the parallel alignment and the filled 
symbols correspond to the antiparallel alignment. In both cases we calculated the superconducting 
in the middle of the S layer. 



on an impurity, tunnels in a ferromagnetic electrode, tunnels back in the superconductor and 
scatters on the same impurity. In the case of strong ferromagnets we find that that a small 
disorder tends to increase the value of T^f , which means a reduction of pair-breaking effects 
(analogous to a reduction of t). This is absent in the antiparallel case where disorder reduces 
the value of T^^ (see Fig. |21-(a)). For larger values of disorder we have T^^ ~ . For weak 
ferromagnets we carried out a simulation at a temperature close to the minimum on Fig.^(b). 
We find that a small disorder tends to increase the value of the superconducting gap in the 
parallel and antiparallel alignments (see Fig.[21-(b)). The increase of the superconducting gap 
is larger in the parallel alignment so that there exists a value of disorder above which the 
superconducting gap in the parallel alignment is larger than the superconducting gap in the 
antiparallel alignment (see Fig.[21-(b)). For larger values of disorder we have Ap(r) ~ Aap(T). 

Effect of a finite thickness in the S layer. - The model of trilayer can be generalized to 
incorporate a finite thickness in the superconductor. This is done by inverting numerically 
the Dyson matrix associated to a chain of L coupled superconducting layers ao : ctL-i- 
The layer and a^+i are coupled by a tunnel amplitude t' . The layer ao (aL_i) is coupled 
to the left (right) to a ferromagnetic layer a (b) by a tunnel amplitude t. 

We have shown on Fig.|31-(a) the variations of the superconducting gap for different values 
of the thickness of the superconducting layer. We have Ap(r) < Aap(T) meaning that pair 
breaking dominates at any temperature for sufficiently large thicknesses. In the case of weak 
ferromagnets with hex of order Aq we find also Ap(r) < Aap(T) for t ^ h^^ (see Fig.|31-(b)). 
We calculated the temperature dependence of the superconducting gap for weak ferromagnets 
with Ls — 2 layers in the superconductor and Lp ^ 2 layers in the ferromagnets and found 
a reentrant behavior. By contrast we found no reentrant behavior for (Ls^Lp) = (2, 1). A 
systematic study of reentrance as a function of Ls and Lp will be presented elsewhere. 
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Conclusions. - To summarize we found that two regimes of parameters are relevant 
from the point of view of possible experiments: strong ferromagnets with tunnel interfaces 
such that Aq <^ t <^ h^^ and weak ferromagnets with tunnel interfaces such that t /lox < Aq 
for F/S/F trilayers with atomic thickness, or /icx — Ao ^ t with a finite thickness in the S 
layer. Quantitative predictions can be obtained from the numerical determination of the 
superconducting gap at finite temperature. For the F/S/F trilayer with atomic thickness and 
with strong ferromagnets we obtain Ap(0) ~ Aap(O) but > T^^ . With weak ferromagnets 
we obtain a reentrant superconducting gap in the F/S bilayer with atomic thickness. We 
discussed the effect of disorder within the simplest approximation in which we neglect vertex 
corrections. In the ferromagnetic alignment and close to the breakdown of superconductivity 
we find that disorder increases in the case of strong ferromagnets, and Ap in the case 
of weak ferromagnets. In the presence of a finite thickness in the superconducting layer pair 
breaking eS'ects dominate and we find Ap (T) < Aap {T) . 

* * * 

The authors acknowledge fruitful discussions with A. Buzdin. 
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